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Abstract 

A slow-rolling scalar field (Quintessence) has been proposed as the origin of 
accelerating universe at present. We discuss some features of quintessence, 
particularly those related with the flatness of its potential. We distinguish 
two types of quintessences, a metrical quintessence which couples to the par- 
ticle physics sector only through the mixing with the spacetime metric and a 
non-metrical one without such property. It is stressed that these two types of 
quintessences have quite different features in regard to the fine tuning of pa- 
rameters required for the flat potential, as well as different implications for the 
quintessence-mediated long range force. The fine tuning argument strongly 
suggests that non-metrical quintessence should correspond to the Goldstone 
boson of an almost exact nonlinear global symmetry whose explicit break- 
ing by quantum gravity effects is highly suppressed. It is briefly discussed 
whether the radion field in higher dimensional theories can be a candidate for 
metrical quintessence. We finally discuss the possibility that string/M theory 
axion corresponds to a Goldstone-type quintessence, and argue that a certain 
combination of the heterotic M theory or Type I string theory axions can be 
a good candidate for quintessence if some conditions on the moduli dynamics 
are satisfied. 
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I. INTRODUCTION 



Recent measurements of the luminosity-red shift relation for Type la supernovae suggest 
that the Universe is accelerating and thus a large fraction of the energy density has negative 
pressure 0. In fact, even before these supernovae data, there have been some indications 
that the Universe is dominated by an yet unknown form of smooth dark energy density 
0. The Big-Bang nucleosynthesis implies the baryon mass density Qp ~ 0.05, while the 
measurements of cluster masses suggest the (clumped) dark matter mass density &dm ~ 0.3. 
(Here all mass densities are expressed in terms of their ratios to the critical mass density 
p c = 3MpH 2 .) On the other hand, the anisotropy of the cosmic microwave background 
shows the first peak of its angular power spectra at I ~ 200, implying that the total energy 
density ^tot ~ 1 which is consistent with the prediction of inflation. This suggests the 
existence of a smooth (unclumped) dark energy density with Qx ~ 0.7, and the supernovae 
data provide a strong support for such dark energy density with further information that it 
has negative pressure density px ~ —0.6Qx- 

If one accepts the accelerating universe as a reality, the utmost question is what is the 
nature of the dark energy with negative pressure. Of course, a nonvanishing cosmological 
constant, i.e. a non-evolving vacuum energy density, would be the simplest form of dark 
energy density with negative pressure. However there has been some prejudice that the 
correct, if exists any, solution to the cosmological constant problem will lead to an exactly 
vanishing vacuum energy density. Recently quintessence in the form of slow-rolling scalar 
field has been proposed as an alternative form of dark energy density with negative pressure 
The true minimum of the quintessence potential is presumed to vanish. However the 
present value of Q is displaced from the true minimum, providing 

V Q ~ HlM 2 P ~ (3 x 10~ 3 eV) 4 , (1) 

and a negative pressure with the equation of state: 

u = ^ = iS-^ < -0.6, (2) 

pq \Q 2 + v Q ~ ' y ' 

where M P = l/y/8irG N = 2.4 x 10 18 GeV denotes the reduced Planck scale and H ~ 10~ 32 
eV is the Hubble expansion rate at present. In fact, when combined with the equation of 
motion for Q, the negative pressure value requires the following slow roll condition: 

(3) 

This condition indicates that if Vq is well approximated by the terms with nontrivial Q- 
dependence, not by a Q-independent constant term, the typical range of Q is of the order 
of Mp or at least not far below Mp. 

The basic idea of quintessence is to assume a scalar field whose potential energy Vq is 
within the order of H^Mp over the Planck scale range of Q. Then an utmost question is how 
such an extremely flat Vq can arise from realistic particle physics models In section II, 
we discuss this issue together with some generic features of quintessence. We distinguish two 
types of quintessences, a Brans-Dicke type metrical quintessence which couples to the particle 
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physics sector only through the mixing with the conformal factor of the spacetime metric 
and a non-metrical one without such property. These two types of quintessences have quite 
different features in regard to the fine tuning of parameters required for the flat potential, as 
well as different implications for the quintessence-mediated long range force. We argue that 
non-metrical quintessence should corrrespond to the Goldstone boson of an almost exact 
nonlinear global symmetry whose explicit breaking by quantum gravity effects is highly 
suppressed. Unless a Goldstone-type, non-metrical quintessence encounters a severe fine 
tuning problem which is much worse than the conventional cosmological constant problem. 
To make this point more clear, we briefly discuss the difficulties of a non-Goldstone type non- 
metrical quintessence having an inverse-power potential Vq = m A+a /Q ot . We also discuss 
whether the radion field in theories with (large) extra dimensions || can be a candidate 
for metrical quintessence [JTU| . In section III, we note that a certain combination of axions 
in heterotic M-theory |TT| or Type I string theory []l2l , |l"3f can be a plausible candidate for 
non-metrical quintessence M if its modulus partner is stabilized at a large value through the 
Kahler stabilization mechanism 
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Although this quintessence axion could avoid the fine 
tunings of parameters associated with its flat potential, it still suffers from the fine tuning of 
cosmological initial conditions. In section IV, we discuss a late time inflation scenario based 
on the modular and CP invariance of the effective potential, which would solve this initial 
condition problem. 



II. SOME GENERIC FEATURES OF QUINTESSENCE 



The supernovae data indicate that the Universe has been accelerating since the redshift 
factor z = 0(1), and thus over the time scale of t ~ l/H ~ 10 10 yrs [|T|. To have a negative 
pressure pq = Q 2 /2 — Vq, one needs 



Q<JV Q ~ H M P 



(4) 



over the accelerating period. In order for the above condition to be satisfied, one needs also 

Q < H%M P (5) 

over the same accelerating period. Note that we have ignored the coefficients of order one, 
and also Q and Q can be regarded as appropriately averaged quantities. At any rate, when 
combined with the equation of motion 



Q + 3HQ + 



9Vq 
dQ 



0. 



(6) 



the Eqs.(|) and (§) lead to 
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(3 x 10- 3 eV)' 
M~ P 



(7) 



which has to be satisfied over the range of Q from the value at z ~ 1 to the present value. 

To see the implications of the slow-roll condition (0), let us assume that Vq is well 
approximated by 
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V Q = V + \m%Q\ (8) 

where V is a Q-independent constant. If Vq is dominated by the constant V during the 
period of interest, Eq.(^) leads to 

\Q/M P \ < (H /m Q ) 2 ~ (1(T 32 eV/m Q ) 2 (9) 

and thus 

\m 2 Q Q 2 /V \ < Hi/ml ~ (1(T 32 eV/m Q ) 2 . (10) 

This shows that if vtlq ^> ifo ~ 10 -32 eV, we are in the limit of non-evolving cosmological 
constant in which Q is frozen at the minimum of Vq. (Here m 2 Q is assumed positive.) 

However if the Q-dependent part i71qQ 2 /2 dominates or at least non- negligible compared 
to the constant Vq, Eq.(|7|) leads to 

\Q\>M P , m Q <l(r 32 eV, (11) 

and so Vq must be extremely flat over the Planck scale range of Q. This is true even when 
mqQ 2 /2 is replaced by other forms of potential, e.g. the inverse power potential m A+a /Q a 
or the exponential potential m i e~ < ^^ v . Thus if it has any feature distinguished from the non- 
evolving cosmological constant, the quintessence model involves an extremely flat potential, 
Vq ~ H^Mp ~ (3 x 10~ 3 eV) 4 over the Planck scale range of Q, and so an almost massless 
scalar boson with mass itlq ~ H ~ 1CT 32 eV in rough order of magnitude estimate. 

As is well known, there are rather strong observational limits on the couplings (and 
also their evolutions) between an almost massless scalar and normal matters |T5[] . For a 
discussion of these constraints and also of the fine tunings required for the flat potential, 
it is convenient to distinguish two types of quintessences: a metrical quintessence which 
couples to the particle physics sector only through the mixing with the spacetime metric, 
and a non-metrical one without such feature. Note that metrical quintessence is quite similar 
to the Brans-Dicke type scalar field in scalar-tensor gravity models. As we will see, these 
two types of quintessences have quite different features in the fine tuning of parameters for 
the flat potential, as well as its implications for the quintessence-mediated long range force. 

Let us first consider non-metrical quintessence Q whose couplings to matter are not 
related with those of the spacetime metric. The interaction lagrangian can include 

Ant = ^-(&* V" + ^G%G^ + C q m q qql (12) 

M P q 

where F^ u and Gf^, denote the electromagnetic and gluon field strengths, respectively, and 
q stands for the light quarks with mass m q . Note that £/s represent the ratios of the 
quintessence couplings to the gravitational coupling. Since the operator combination which 
couples to Q does not match to the energy momentum tensor even in the non-relativistic 
limit, they violate the equivalence principle for non-relativistic laboratory bodies. More 
explicitly, the relevant operator combination can be written as 

faG^G^ + Y.^ m q qq = ^CPTH + - J^j^ + l m )}m q qq, (13) 
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where is the trace of the QCD energy momentum tensor with the QCD beta function 
(3(g) and the mass anomalous dimension j m . Then the second term in the r.h.s. gives rise 
to a composition-dependent acceleration of the i-th test body, yielding 



where £ 3 = £ g — 2g£ G 2(l + r y m )/^(g) and the subscript "s" denotes the source of the long 
range force. Under the assumption that there is no sizable cancellation between £ q and 
2fi , £G 2 (l+7m)//3(<?), which is the case for non-metrical quintessence, the current observational 



limit on non- universal acceleration [15 



|Ao/o| < 1(T 12 (15) 



leads to 16 



£ G2 < 0(1O- 5 ), Z q <O(10- 5 ). (16) 

for the operator coefficients £q2 and £ q renormalized at energy scale ~ 1 GeV. 

If t;p2 is nonzero at low energies, say at 0(1) eV, a time- varying Q leads to a time- varying 
fine structure constant which is constrained [01 as 



aem/a em < 10- 16 yr _1 - (17) 

Obviously this gives 

< 3 x l(T 6 Mp# , (18) 

a rather strong observational limit on the combination ^f 2 Q- 

The direct observational limits (0) already suggest that the non-derivative couplings 
between non-metrical Q and normal matter are highly suppressed, even compared to the 
gravitational strength. In fact, although indirect, much stronger limits on generic non- 
derivative couplings can be obtained from the fine tuning argument associated with the 
extremely flat potential. The quintessence couplings of Eq. ([12]) mean that the electromag- 
netic and QCD fine structure constants, and the light quark masses are all Q-dependent. 
Then the QCD and electroweak physics contribute to Vq, yielding the following form of the 
quintessence potential in rough order of magnitude estimates 

Vq = Kcd + m q A QCD + jl^M^ ln(M 2 w /^) + ... (19) 

where the ellipsis stands for other types of contributions to Vq. Here the QCD scale Aqcd oc 
e -2Tr/ba QC D{Q) anc [ the W-mass Mw oc a em (Q) are some functions of Q whose Q-dependence 
would be determined by the Q couplings to the gluon and electroweak field strengths. 

Suppose we tuned the parameters of the theory to make Vq ~ (3 x lCT 3 eV) 4 for a 
particular value of Q, say for Q = Q . Then for different value of Q, e.g. Q = Q + AQ, 
the three terms in Vq of Eq. fll9D would lead to a potential energy difference 
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S 2 V q ~ e ? m 9 A| CD ^ ~ io 42 e^ 2 M|,^ 



6 3 V Q ~ 10- 2 ^M^ ~ io 52 ^ 2 ^ 2 M|^. (20) 

These estimates of the variation of Vq show that even after the fine tunining for Vq ~ H^Mp 
at Q — Qq, one still needs additional fine tunings in order to keep Vq within the order of 
H^Mp over the Planck scale range of Q. Generically, one needs different fine tunings for 
different values of Q, and then the theory would suffer from the infinite numbers of fine 
tuning problems. 

There are essentially two ways to avoid this difficulty The first one is to assume a non- 
linear continuous global symmetry (U(1)q) under which Q is shifted by a constant. Such 
U(1)q would assure that all non-derivative couplings of Q to the particle physics sector are 
suppressed, so that the particle physics sector contribution to Vq can be small enough for 
the wide range of Q. The second possibility is that although the non- derivative couplings 
between Q and normal matter are sizable, they are all precisely correlated to each other in 
such a way that their contributions to Vq are cancelled independently of the values of Q. 
It is likely that the only sensible model in this direction is the metrical quintessence whose 
couplings to matter are induced only through the mixing with the spacetime metric. This 
leads us to conclude that, in order to avoid the infinite numbers of fine tuning problems, 
non-metrical quintessence should correspond to the Goldstone boson of a nonlinear global 
symmetry under which 

U Q (1) : Q^Q + constant. (21) 

Still the potential difficulty for Goldstone-type quintessence is that it requires a nonzero 
but extremely tiny breaking of U(1)q. Note that generically a global symmetry can not be 
arbitrarily good, particularly when quantum gravity effects are included. So the questions 
are (1) how U(1)q could avoid a potentially large breaking from quantum gravity effects and 
(2) what is the origin of the tiny U(1)q breaking which would provide small but nonvanishing 
Vq ~ H^Mp. To quantify this issue, let us first estimate the allowed size of U(1)q breaking. 
Again assuming that U(1)q breaking couplings of Eq . (|T2"D are not correlated to each other, 
the variations of Vq in Eq. ( p0| ) suggest that 

0,2 <Hn 52 , £ G2 <10- 43 , ^ < lO" 42 , (22) 

in order for Vq remain to be within the order of H^Mp over the Planck scale range of Q. 
This already shows that [7(1 ^-breaking couplings should be extremely suppressed, even 
compared to the gravitational coupling strength. 

In fact, the variation ( p0|) of Vq includes only the contributions from the particle physics 
degrees of freedom with momenta p ^ My/- If one includes the contributions from higher 
momentum modes, the typical variation of Vq takes much larger value. As an illustrative 
example, let us consider generic supergravity models with cutoff scale A ~ Mp. The Kahler 
potential K, superpotential W, and gauge kinetic functions f a can be expanded in powers 
of the visible sector matter fields $*: 
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K = K {Z, Z*) + Zij(Z, Z*)&$* j + Hij(Z, Z*)&& + X ijk (Z, Z*)¥^>^* k + 



W = W (Z) + Y ijk (Z)$ l &$ k + T ijkl (Z)<$> l &<& k <& 1 + 
f a = f a0 (Z) + f aij (Z)&& + ..., 



(23) 



where Z denotes the chiral multiplet including the quintessence degree of freedom Q. Any Q- 
dependence of the coefficient functions in the expansion gives rise to non-derivative couplings 
between Q and $ l , so reflects the explicit breaking of U(1)q. The typical sizes of the 
variation of Vq induced by each of the Q-dependent coefficient functions are estimated in 
||, taking into account various types of quantum corrections. Again if any of the non- 
derivative couplings is sizable, we can not keep Vq small enough over the Planck scale range 
of Q by tuning finite number of parameters. This consideration leads to the following strong 
limits on the U(1)q breaking couplings (in the unit with Mp = 1): 



\n(K ) < 10 



JL Z . .< io- 81 — 
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d_ 
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-74 



Yi 



ijk 



I 2 < 



10 



-79 



df 



at J 



< 10 



-67 



dz t 

lm(f QCD ) < 10^ 42 , 



(24) 



where we have assumed an weak scale gravitino mass and the last bound is on the coupling 
of Q to the QCD anomaly G^G^. In fact, even the flat potential argument does not give 
a meaningful limit on the coupling to the electroweak anomaly. So it is in principle possible 
that £pp is sizable for the coupling 



Q 



(25) 



Unless \£ FF \ 10 -2 , such coupling may lead to an observable rotation of the polarization 
axis of the radiation from remote sources [16 , 17] . However in view of the gauge coupling 
unification, it is somewhat unlikely that Q has a sizable coupling to FF, while its coupling 



to GG is so small as given by the last bound of Eq. fl24j) . 

The above limits on the U(1)q breaking couplings suggest that U(1)q must be an almost 
exact global symmetry. It is highly nontrivial to have such an almost exact global symmetry 
in any realistic models for particle physics once the effects of Planck scale physics are included 
through arbitrary higher dimensional operators in the effective action. In the next section, we 
will explore the possibility of such global symmetry in the framework of string/M theory and 
point out that a certain combination of the heterotic M theory or Type I string theory axions 
can be a plausible candidate for quintessence if some conditions on the moduli dynamics 
are satisfied. Here we stress that U(1)q has nothing to do with a small Vq at its minimum. 
It just ensures that Vq is small over the Planck scale range of Q once it were small at 
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its minimum. Without U(1)q, non-metrical quintessence model would suffer from infinite 
numbers of fine tuning problems to keep Vq to be small over the wide range of Q. At any 
rate, the limits ( p2[ ) which are based on the fine tuning argument suggest that it is hopeless 
to see the long range force mediated by non-metrical quintessence. 

Let us now discuss the Brans-Dicke type metrical quintessence which couples to matter 
only through the mixing with the conformal factor of the spacetime metric. We first consider 
the observational limits on its couplings. The dynamics of metrical quintessence would be 
described by the following effective action: 

S eff = J d A x [v=£ (-\m 2 p k{ 9 ) + ^ u d,Qd u Q - Vo(Q)) + \pg £m(g^, $)] (26) 

where 1Z is the Ricci scalar density and L m denotes the lagrangian density of generic matter 
fields $ (including the conventional particle physics sector) under the background spacetime 
metric given by 

~g, v = A 2 (Q/M P )g^. (27) 
Then the couplings between Q and $ are given by 

d ( r~~r \ <- Q 



c '-= Q sQ{^H Q , Q r^ < 28 > 



where 



'-'-mi,. 

and 

T£ = ^G a ^G% + ^-F^F^ + ]T(1 + lm )mMi + - (30) 

denotes the trace of the energy momentum tensor of the matter fields $ = G^ v , F^,^. 
Since it respects the equivalence principle in non-relativistic limit, the above quintessence 
couplings lead to a deviation from the Einstein's general relativity only through the relativis- 
ts corrections. Currently performed gravitational experiments in the solar system provide 
upper bounds on such deviation |15| , implying 

£<4.5xl0- 2 . (31) 



One might think that the strong observational limit fllTD on the time-varying fine structure 
constant can give a constraint on the evolution of metrical quintessence. However existing 
limits are derived from the electromagnetic processes with momentum transfer much smaller 
than the electron mass. For such low momentum transfer, we have /3(e) = and so a time- 
varying metrical quintessence does not lead to a time-varying electromagnetic fine structure 
constant. 

Let us now turn to the flatness of the metrical quintessence potential. Since the couplings 
between Q and the particle physics sector are determined by the mixing with the spacetime 
metric, the whole contribution to Vq from the particle physics sector is simply given by 
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AS eff = J d 4 Xy/^AV Q = J d'x^gKt = J d 4 v ^A 4 (Q/M P )A 4 mJ (32) 

where is a Q -independent energy density which is a generic function of all particle physics 
parameters in C m . Note that the general covariance w.r.t guarantees the above form 
of AVq which is obtained after integrating out all particle physics sector fields. The total 
quintessence potential energy is then given by 

V Q = V + AV q = V (Q/M P ) + A\Q/M P )A i m . (33) 

An interesting feature of the metrical quintessence potential is that once the particle 
physics parameters are adjusted to yield A^ < HqM p , the entire contribution to Vq from 
the particle physics sector remains to be within the order of H^Mp over the Planck scale 
range of Q. Of course one still needs to adjust the theory to keep Vo to be small enough. 
The typical mass scales of the physics responsible for Vo can be much lower than those of the 
conventional particle physics and/or the physics generating Vq may enjoy extra symmetry 
which would assure that Vo is small enough. In this case, one may be able to make the total 
Vq = Vq + AVq flat enough without a severe fine tuning problem other than the minimal fine 
tuning for small A^. Note that if the (non-derivative) coupling between metrical quintessence 
and normal matter takes a value not far below the current observational limit ([H]), e.g. 
£ = C(10 -2 ), it can lead to observable effects in future gravity experiments. 

An interesting candidate for metrical quintessence would be the radion field in theories 
with large extra dimensions [0,0]- However at least in the simplest framework, the radion 
has too strong coupling with matter to be a quintessence. For example, if one considers a 
higher dimensional theory with n large flat dimensions with common radius, the canonical 
radion field Q in the Einstein frame of the four dimensional metric appears in the 
(4 + n)-dimensional metric: 

ds\ i+n) = g£ + N n) dx M dx N = et Q / Mp 9llv dx»dx v + e^^dyW (34) 

where 

£ = yV^ + 2). (35) 

If the whole particle physics sector lives on a three brane, the radion couples to the particle 
physics sector only through the induced metric on the brane: 

g $m) = e ^Q/M Pg ^ (36) 

Obvioulsy, the value of £ in this framework does not satisfy the observational limit (PIP , 
and thus the radion Q can not be a quintessence. (Of course, once one makes the radion 
massive by suitable stabilization mechanism fl8 |, the limit ([31] ) can be avoided.) Although a 



Brans-Dicke type quintessence can not arise from the simplest version, it may be possible to 
arrange the higher dimensional theory and the compactification scheme to yield an almost 
massless radion field with £ < O(10~ 2 ). 

Let us close this section with a brief discussion of the difficulties of specific non-metrical 
quintessence which is not Goldstone type, but has an inverse power law potential: Vq ~ 
m 4+a /Q a . Perhaps the most attractive feature of this model would be the tracker behavior 



9 



that pq ~ Pmatter at present for wide range of initial conditions of Q and Q. This feature 
allows us to avoid the cosmic coincidence puzzle asking why now pq is comparable to Pmatter 
although their cosmological evolutions are generically different H||. 

It has been pointed out 0] that the inverse power potential can be generated by non- 
perturbative dynamics in the underlying supersymmetric theory pl| . The model includes 
a hidden supersymmetric QCD sector with SU(N C ) gauge group and Nf flavors of hidden 
quarks and antiquarks $ = ($;, It is then assumed that the tree level superpotential is 
independent of $. The ^-dependent part of the exact superpotential is then given by 

/ \3Nc-Nf \1/(N C -N f ) 

w *- = (sw) • (37) 

where A denotes the dynamical scale of the hidden SU(N C ) gauge group. If Wdyn is a 
dominant term in the full superpotential, the effective potential of the SU(Nf) invariant 
D-flat direction $j = $^ = Q is given by 

V Q = FA A+a /Q a , (38) 

where a = 2(N C + Nf)/(N C — Nf) and T is a function of Q/Mp which is essentially of 
order unity. The explicit form of T depends on the Kahler potential of $. For A ~ 
M P (H /Mp) 2/{4+a) ~ 10~ 120/(4+a) M P , this potential is within the range of H$Mp over the 
Planck scale range of Q, so can be a candidate for quintessence potential. 

However the problem is that Wd yn ~ H Mp <^ m 3 / 2 Mp for Q ~ Mp, so it can not be 
a dominant term in the superpotential in any realsitic model of supersymmetry breaking. 
For Q ~ M P , supergravity effects become important, and so we have to consider Vq in the 
supergravity framework in which the scalar potential is schematically given by 

Vsugra = \F\ + D , (39) 

where F and D denote the auxiliary F and D components of supersymmetry breaking fields, 
and W is the total superpotential. For realistic phenomenology, we need |F| > (and/or 
D > M$y) for the weak scale mass M w . Obviously we then need W ~ M P \F\ ~ m 3 / 2 Mp 
in order to have the total vacuum energy density much smaller than M^. Since Wd yn of 
Eq.(^7)) is just of the order of H Mp, the total superpotential must include an additional 
term W ~ m 3 / 2 Mp. Once W is introduced, which is a necessity, then the flatness of Vq is 
totally lost. 

To see this explicitly, let us consider the case with 

W = W + W dm + e 1 { -^T + ..., 



Mp 



K = + + e^^Sr + (40) 



Mp 



where n and k are some integers which can be arbitrarily large, and the ellipses stand for 
other possible terms. We then have (again schematically) 
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j^4+a Q 2 ^ Q 2n Q^ n ~ 2 

V Q = -Q^T + m l/2Q 2 + e 2 m l/2jj2k^2 + ^3/2^2^3 + ^M^=S + -» ( 41 ) 



Obviously, for any realistic value of 7713/2, even after tuning the theory to have Vq < H^Mp 
at its minimum, one can not keep this Vq within the range of H^Mp over the Planck scale 
range of Q. 



III. HETEROTIC M OR TYPE I STRING AXION AS QUINTESSENCE 

In the previous section, we argued that non-metrical quintessence should corrrespond 
to the Goldstone boson of an almost exact nonlinear global symmetry in order to avoid 
infinite number of fine tuning problems associated with the flat potential. It is also noted 
that explicit breaking of this global symmetry by quantum gravity effects must be extremely 
tiny. In this section, we discuss the possibility that string and/or M theory axions correspond 
to such Goldstone-type quintessence. There are three different classes of string/M theory 
vacua which are particularly interesting in phenomenological sense: perturbative heterotic 
string vacua [I3|, heterotic M theory vacua [IT], and Type I string vacua with D-brane 



configurations [0. It is well known that axions in perturbative heterotic string vacua 
receives a large potential energy associated with the string world sheet instanton effects 
T9fl , so here we concentrate on heterotic M theory and Type I string vacua. We will see 



that certain combination of the heterotic M theory or Type I string theory axions can be a 
plausible candidate for quintessence if some conditions on the moduli dynamics are satisfied. 

Let us first consider the heterotic M-theory on a 11-dimensional manifold with boundary 
which is invariant under the ^-parity ||11||: 



C^-C, z u ->-z u , (42) 

where C = CABcdx A dx B dx c is the 3-form field in the 11-dimensional supergravity. When 
compactified to 4-dimensions, axions arise as the massless modes of C^n and C mn \\. 

^ upa d [u C pa]11 = d^ S , C m nll = E Vi(^)<n, (43) 

i 

where uo l (i = 1 to hi t \) form the basis of the integer (1,1) cohomology of the internal 6- 
manifold and fi, v are tangent to the noncompact 4-dimensional spacetime. In 4-dimensional 
effective supergravity, these axions appear as the pseudo-scalar components of chiral multi- 
plets: 

S =(Att)- 2 / 3 k^ 3 V + t7] S} 

Ti = (47r)- 1 / 3 /t- 2 / 3 f uA dx 11 + i-m, (44) 
Jd 

where k 2 denotes the 11-dimensional gravitational coupling, V is the volume of the internal 
6-manifold with the Kahler two form a>, and the integral is over the 11-th segment and also 
over the 2-cycle Cj dual to uj 1 . Here the axion components are normalized by the discrete 
Peccei-Quinn (PQ) symmetries: 

Im(5) -»■ Im(S') + 1, Im(Ti) -> hn(ty + 1, (45) 
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which are the parts of discrete modular symmetries. 

Holomorphy and the discrete PQ symmetries imply that in the large Re(S') and Re(Tj) 
limits the gauge kinetic functions can be written as 

4vr/ a = k a S + l*Ti + 0(e-^ s , e" 2 ^), (46) 

i 

where k a and l ai are model-dependent quantized real constants and the exponentially sup- 
pressed terms are possibly due to the membrane or 5-brane instantons. For a wide class of 



compactifed heterotic M-theory, we have [22,23 



i 

4tt/ e , =S-J2hT l + 0(e- 2 * s , e" 2 ^), (47) 

i 

where ZjT, corresponds to the one-loop threshold correction in perturbative heterotic string 
terminology with the quantized coefficients determined by the instanton numbers on the 
hidden wall and also the orbifold twists. 

Let Q denote a linear combination of Im(Tj) orthogonal to the combination Ya ^lm(Tj) 
and define the associated nonlinear global symmetry: 

U(1) Q : Q^Q + constant. (48) 

Note that such Q exists always as long as hn > 1. In this regard, models of particular 
interest are the recently discovered threshold-free models with /j = |23| for which any of 
Im(Tj) can be identified as Q. At any rate, Q can be a candidate for quintessence if U(1)q 
breaking couplings in the effective supergravity are all suppressed as in Eq. 
It is easy to see that for Q defined as above 

®fa /nf .-1-kT 



dQ 



e>(e- 27rJ ), (49) 



over the entire range of Q. Here we use the unit with Mp = 1 and the internal 6-manifold 
is assumed to be isotropic, so Re(Tj) ~ Re(T) for all Tj. Holomorphy and discrete PQ 
symmetries imply also 



dQ 



0{e~ 27TT ) (50) 



for the Yukawa couplings in the superpotential. Similar estimate applies also to the 
non-renormalizable holomorphic couplings in the superpotential. A non-perturbative super- 
potential Wq may be induced by these holomorphic gauge and superpotential couplings, e.g. 
by gaugino condensation. The U(1)q breaking in W is estimated to be 

aw dw df a dW dY m „,_„„,, . 

We then have 
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(52) 



for the full superpotential W. 

The above estimates of U(1)q breaking couplings in /„ and W are made by simple macro- 
scopic argument based on supersymmetry and discrete PQ symmetries. However one can 
easily identify its microscopic origin by noting that 27rRe(Ti) = (47r) -1 / 3 /t~ 2 / 3 j c . u> A dx 11 
corresponds to the Euclidean action of the membrane instanton wrapping C, and stretched 
along the 11-th segment f2"2]| . When extrapolated to the perturbative heterotic string vacua, 



such membrane instanton corresponds to the heterotic string worldsheet instanton wrap- 
ping the same 2-cycle ||19|| . Explicit computations then show that the Kahler potential 



and/or the gauge kinetic functions are indeed corrected by worldsheet instantons, yielding 
5K = 0(e~ 2nT ) and 5f a = 0(e~ 2nT ) |24]]. These corrections can be smoothly extrapo- 
lated back to the heterotic M-theory vacua |26| and identified as the corrections induced 
by stretched membrane instanton. If a nonperturbative superpotential Wq is generated by 
gaugino condensation, 5f a = 0(e~ 2lTT ) leaves its trace in 6W = O(e~ 2lTT W ). We thus 
conclude that U{1)q breaking terms of C(e~ 2?rT ) are indeed induced in K, W and f a for 
generic heterotic M theory vacua. 

The holomorphy and discrete PQ symmetries ( [45]) ensure that U(1)q breaking terms in 
f a and W, whatever their microscopic origin is, are all suppressed by e~ 2nT , Also the U(1)q 
breaking terms in K induced by stretched membrane instanton are suppressed by e~ 2nT , 
However there may be unsuppressed U(1)q breaking term in non-holomorphic K, which 
would arise from yet unknown microscopic origin. Although not definite, it is unlikely to 
have such unsuppressed correction. U{1)q is a linear combination of U{l)i : <5Tj = ici (cj = 
real constant) which originate from the local transformation of the 11-dimensional three form 
field: 5C = w'A dx 11 with q = / 5C = J uj % A dx 11 . This indicates that upon ignoring the 
effects of boundary degrees of freedom, the zero momentum mode of Im(Tj) couples only 
to stretched membrane instantons wrapping the 2-cycle Ci and are stretched between the 
boundaries. Since we chose the combination to avoid the breaking by Yang-Mills instantons 
on the boundary, U(1)q appears to be broken only by stretched membrane instantons whose 
effects are suppressed by e~ 27rT . It is thus expected that 

% = ( e-^) (53) 

as in the case of f a and W. 

Obviouslly Eqs. (f||), (|52D , and (]53D show that U(1)q becomes an almost exact global 
symmetry in the limit Re(T) ^> 1. From the supergravity potential 



V = e 

' sugra ' 



K 



K IJ DjW(DjWy - 3\W\ 2 (54) 



and also the standard order of magnitude relations m 3 / 2 ~ W/Mp ~ DiW/Mp, one easily 
finds that the axion potential is given by 

V Q ~ e- 2wReiT) ml /2 Mp cos[27rIm(T)]. (55) 

This axion potential can be identified as the quintessence potential Vq ~ (3 x 10~ 3 eV) 4 if 
the modulus vacuum value is given by 
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Re(T) ~ i- ln(m^ /2 M|,/F Q ) ~ 32, (56) 

where TO3/2 ~ 10 2 GeV are used for numerical estimate. 

In the above discussion, gauge kinetic functions are assumed to be given by fl4"7|). In such 
cases, Q must be a linear combination of Im(Tj) to avoid the couplings to the QCD and/or 
the hidden gauge anomaly. Gauge kinetic functions can be generalized to the form of fljSp 
with non-universal k a , and then Q can include the Im(S) component. At any rate, as long 
as the quantized coefficients k a and l a i are all of order unity, the domain of moduli space 
allowing a quintessence axion in heterotic M theory has 

Re(S) = O (^—) , Re(T) = O (^—) ■ (57) 

On this domain, we also have 

Mp 

K 2/3 



; [Re(5)Re(T)]- 1 / 2 = 0{a GUT ), 

Re(S)/[Re(T)] 3 = 0(a 2 GUT ), (58) 



implying that this domain gives the unification scale Mqut close to the phenomenologically 
favored value 3 x 10 16 GeV, and also the length np of the 11-th segment about one order of 
magnitude bigger than the 11-dimensional Planck length /t 2//9 . 

Let us now turn to Type I string axions. The Type I axions (again normalized by the 
discrete PQ symmetries of fl4"5"D) correspond to the massless modes of the R-R two form 



fields B^ u and B mn : 



and they form 4-dimensional chiral multiplets together with the string dilaton e D and the 
internal space volume V: 

S = (2n)- 6 a'- 3 e- D V + i Vs , 

Ti = (27r)- 2 a'- l e- D [ uj + i^. (60) 
Jd 

Here we include D9 and D5 branes in the vaccum configuration, and consider the 4- 
dimensional gauge couplings a$ and defined on D9 branes wrapping the internal 6- 
manifold and D5 branes wrapping the 2-cycles Cj, respectively ]l2j. Again holomorphy and 



discrete PQ symmetries imply that the corresponding gauge kinetic functions can be written 
as ([46]). A simple leading order calculation gives a 9 = l/Re(S) and a 5i = l/Re(Tj), and so 

m 



Anf 9 = S, ^hi = Ti. (61) 

However this leading order result can receive perturbative and/or non-perturbative correc- 
tions. Generic perturbative corrections can be expanded in powers of the string coupling 
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e D , the string inverse tension a', and also the inverse tension e D a' k of -C^fc-i branes [13 . 
Generically they scale as 

e nD a' m oc [Re(^)]^[Re(T)]^. (62) 



Combined with the leading order result (pl|) and also the general form of gauge kinetic func- 
tion (f46|) dictated by holomorphy and discrete PQ symmetries, this scaling behavior implies 
that /g can receive a Tj-dependent correction at order a 12 , while there is no perturbative 
correction to fa, and so 

4tt/ 9 = S + UTi + 0(e- 2nS , e~ 2 ^), 

47r/ 5 ,=T, + 0(e~ 2 ^,e~ 27rTl ). (63) 

Similarly to the case of heterotic M-theory, the quintessence axion Q can arise as a linear 
combination of Im(S') and Im(Tj), however its explicit form depends on how the various 
gauge couplings are embedded in the model. Note that for the QCD and even stronger 
hidden sector gauge couplings, both Re(/ a ) and Im(/ ) are required to be Q-independent 
as in fl24]), while for the weaker gauge interactions Im(/ a ) are allowed to have a sizable In- 
dependence. Here are some possibilities. If all of the standard model gauge couplings and 
the stronger hidden sector gauge couplings are embedded in a 5i , lm(S) can be a quintessence 
when Re(S') ~ 32. If those gauge couplings are embedded in a 9 , any of Im(Tj) can be a 
quintessence when Re(Tj) ~ 32. If the vacuum does not include any D5 brane wrapping the 
particular 2-cycle, e.g. the i-th cycle, the corresponding axion Im(Tj) can be a quintessence 
independently of the embedding when Re(Tj) ~ 32. 

For the quintessence component Q defined as above, the Q-dependences of K, W and 
f a are all suppressed by e~ 2wZ (Z = S or T) as in the case of heterotic M-theory. So 
the quintessence axion potential is again given by ([55]) where now T is replaced by Z. 
Microscopic origin of this Type I axion potential can be easily identified also by noting that 
27rRe(S') corresponds to the Euclidean action of D5 brane instanton wrapping the internal 
6-manifold and 27rRe(Tj) is of Dl string instanton wrapping the 2-cycle C«. Quintessence 
axion in Type I string theory requires also that both Re(S') and Re(T) are of 0(1/ o>cut)- 
On this domain of moduli space, we have 

e 2D ~ Re(5)/[Re(T)] 3 = 0(a 2 GUT ), 

a'/V 1/s ~ [Re(T)/Re(5)] 1/2 = 0(1), (64) 

and thus a rather weak string coupling and also a rather strong sigma model coupling. 

So far, we have noted that a certain combination of axions in heterotic M theory or Type 
I string theory can be a candidate for quintessence if its modulus partner has a large vacuum 
value. For the axion component normalized as lm(Z) = Im(Z) + 1, explicit breaking of the 
associated global symmetry (U(1)q : lm(Z) — > lm(Z) + constant) is suppressed by e~ 2nZ . 
More concretely, the effective supergravity model is described by 

K = K(Z + Z*) + 0(e- 2liZ ), 

fa = fa + 0(e- 2 * Z ), 

W = [1 + 0(e~ 27TZ )]W, (65) 
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where f a denote the gauge kinetic functions for the standard model gauge group and also 
the stronger hidden sector gauge group, and f a and W are Z -independent. We then have 

£q = \{d,Qf - m 4 [cos(Q/v Q ) + 1], (66) 

where 

m 4 ~ e- 2 ^ z W 3/2 M 2 P ~ (3 x KT 3 eV) 4 , (67) 

for Re(Z) ~ 32, and the canonical quintessence axion Q and its decay constant vq are given 
by 

Q = M P V2K"lm(Z), v Q = —M P V2K", (68) 

where K" denotes the Kahler metric of Z: 

K" = (69) 

At leading order approximation in perturbation theory, we have K = cln(Z + Z*) for 
a constant c of order unity, and then Re(Z) can not be stabilized at the desired value. 
Obviously W is almost Z-independent for large Re(Z), so can not provide a stabilizing 
potential of Re(Z). Note that if W provides any sizable potential of Re(Z), it means also a 
sizable potential of lm(Z), which is not allowed for Im(Z) to be a quintessence. We thus need 
a rather strong ?7(l)g-preserving nonperturbative effects encoded in the Kahler potential, 
e.g. K np ~ (Z + Z*) k e~ wz+Zm , which would stabilize Re(Z) at Re(Z) ~ 32 @. As we will 



discuss in the next section, such strong correction to K is required also to enlarge K" for 
cosmological reasons. 



IV. INITIAL CONDITION PROBLEM AND A LATE TIME INFLATION 

SOLUTION: 

In order to suppress U(1)q breaking quantum gravity effects, we needed a large vacuum 
value of Re(Z). A large value of Re(Z) then implies K" = 0([Re(Z)]~ 2 ) 1, and so 
vq <C Mp. (See Eq. fl68|).) Since Vq determines the coupling strength of Q, while Mp 
determines the gravitational coupling strength, Q responds to its potential energy more 
sensitively than the expanding universe does. As a result, generic initial values of Q and Q 
give a rapidly rolling Q and thus positive pressure at present. To avoid this, we need a fine 
tuning of initial condition. So the attempt to avoid the fine tunings of parameters for the 
flat potential leads to a new fine tuning problem. In this section, we propose a late time 
inflation scenario based on the modular and CP invariance which would solve this initial 
condition problem. 

Let us first consider what kind of initial conditions we need to have. When applied for 
( |55D and (|6"8"D, the slow roll condition (|7]) leads to 
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|27rIm(Z)| < V2K"/2n = O 



27rRe(Z) j 



|27rIm(Z)| < 2nH /V2K^ = O (|^|) . (70) 

implying that Im(Z) should be at near the top of its potential. The main difficulty of this 
condition is that 



m Q 



d 2 v 



dQ 2 



1/2 ^ V3 



™ « ^Li/ » if , (71) 



and thus a small Im(Z) is unstable against the cosmlogical evolution during the period with 
an expansion rate H ^ tuq. The resulting instability factor is given by e yrnQ ^ H ° where 7 is a 
constant of order unity whose precise value depends upon the initial conditions. A detailed 
numerical study of the cosmological evolution gives 7 ~ 0.5 for wide range of relevant initial 
conditions Thus the slow-roll condition fl70"| ) requires eventually the following initial 

conditions 

\2nlm(Z)\ m < e -W"o „ e -3/^ 

\2irlm{Z) | in < He~^ /Ho ~ He~ 3/VW ' (72) 

for the period with if 3> mg. 

Eq. ([72]) shows that the degree of required fine tuning is quite sensitive to the value of 



the Kahler metric K" . At leading order approximation in string or M theory, we have [13 
K pa —cln(Z+Z*) where c is a constant of order unity, so K" = O(10~ 3 ) for Re(Z) ~ 30. As 
we will see, the degree of fine tuning for this value of K" is too severe to be accommodated, so 
we need a mechanism to enlarge the value of K" . It is expected that K" can be enlarged by 
?7(l)Q-preserving nonperturbative effects which would be responsible for stabilizing Re(Z) 
|T^| . However it is hard to imagine that K" becomes of order unity, which would allow to 
avoid the fine tuning of initial condition without any further mechanism. In the following, 
we discuss a late time inflation scenario which would resolve the fine tuning problem of 
initial conditions for a reasonably enlarged value of K". 



The gauge symmetries of string or M-theory include discrete modular group JT3 under 
which Z and other generic moduli transform as 

Re(Z) ^RelX>' MZ)^lm(Z) + l, 0-0', (73) 
and also CP fE7| under which 



Z — Z*, — 0*. (74) 

Here we take the simplest form of the Z-duality [Z = S or T), i.e. Re(Z) — l/Re(Z) 
with the self-dual value Re(Z) = 1, however our discussion is valid for other forms of the 
Z-duality transformation as long as the self-dual value is of order unity. 
The modular and CP invariance ensure that the invariant points, 
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Re(Z) = 1, lm{Z) = or J, = 0* or 0'*, (75) 



correspond to the stationary points of the effective action [pq] . It is then quite possible that 
the modular invariant point Z = 1 is a (local) minimum of the effective potential during the 
inflationary period if the inflaton is a modular invariant field. This (local) minimum may 
become unstable if the inflaton field takes the present value, making Ke(Z) rolls toward the 
present minimum at Ke(Z) ~ 32 after inflation. At any rate, during the inflationary phase, 
Ke(Z) ~ 1 and then all the moduli masses including that of lm(Z) have the same order of 
magnitude: 

mRe(z) ~ mm(z) ~ m<f, ~ H in f, (76) 

where Hi n f denotes the expansion rate during inflation. Note that the axion mass m\ m ^z) is 
unsuppressed for Re(Z) ~ 1. Since Re(Z) was far away from the present value, it is expected 
that the inflationary potential is at least of 0(my 2 Mp), and so Hi n f is at least of 0(^3/2)- 
To avoid a too large quantum fluctuation during this inflation, we take the minimal value 
of H inf , so 

H inf = 0(m 3/2 ). (77) 

About the location of the minimum in the axion direction, we have just two possibilities 
if CP is not spontaneously broken in the moduli sector. One of the two CP invariant points, 
lm(Z) = and 1/2, is the minimum, while the other is the maximum. Our key assumption 
is that lm(Z) = was the minimum for the inflationary modulus value Ke(Z) = 1, however 
it becomes the maximum for the present modulus value Ke(Z) ~ 32. Note that the coefficient 
of the cosine potential of lm(Z) is a function of Re(Z), and so its sign can be changed when 
Ke(Z) varies from the inflationary value to the present value. 

Given the features of the moduli potential discussed above, during the inflationary period 
all moduli are settled down near at the modular and CP-invariant local minimum with 
Ke(Z) fa 1, lm(Z) w 0, and w 0*. In particular, just after inflation, we have 



|27rIm(Z)U = 0{e~^' 2 ) + 0(H inf /M P ), 

\2nlm(Z)\ in = 0(e- 3N ^ 2 H mf ) + 0(H? nf /M P ), (78) 

where N e denotes the number of efoldings. Here the exponential suppression is due to the 
classical evolution toward the minimum of the inflationary potential, while the second terms 
represent the quantum fluctuations. 

After this late inflation, Re(Z) rolls toward the present minimum at Re(Z) ~ 32. In this 
period, the coefficient of the axion potential changes its sign, and thus lm(Z) = which 
was the minimum at the inflationary phase becomes the maximum of the present axion 
potential. The small value of lm(Z), i.e. ([78]), which was set during the inflation becomes 
unstable if H starts to be smaller than tuq. However if the inflationary values (|78| ) satisfy the 
condition (|72|) , the present value of the quintessence axion satisfies (0), and thus provides 
an accelerating universe at present. This requires a large number of efolding 

> 2JTUQ J_ 
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and also a small quantum fluctuation 



SQ~^< Mpe- 3 ^. (80) 

However the number of efolding can not be arbitrarily large. For a relatively late inflation 
with H in f ~ ?7i3/2, it is extremely difficult to generate the observed density fluctuation 
Sp/p ~ 1CT 5 . It is thus reasonable to assume that density fluctuations were created before 
the late inflation. In this case, the late inflation is required not to destroy the pre-exisiting 
density fluctuations. This consideration gives an upper bound on N e yielding 

N e < 25 ~ 30 (81) 

for a weak scale 7713/2 and the reheat temperature range T r = 10 5 ~ 10~ 2 GeV. Combining 
([79]), (|80|) and fl8~I]) together, we find that our late time inflation scenario can successfully 
generate the required initial condition (|72|) if the Kahler metric is enlarged to be 

K" w 8 x 10 -3 . (82) 

It is not unreasonable to expect that the nonperturbative terms in the Kahler potential 
stabilizing Re(Z) provide the necessary enlargement of the Kahler metric, particularly when 
those terms involve a large power of (Z + Z*), e.g. K np = h(Z + z*) k e~ b ^ z+z * = 0(1) for 
Re(Z) ~ 32 with k > 6, 
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